Single-photon-level quantum memory at room temperature 
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Quantum memories capable of storing single photons are essential building blocks for quantum 
information processing, enabling the storage and transfer of quantum information over long dis- 
tances [TJ [2]. Devices operating at room temperature can be deployed on a large scale and integrated 
into existing photonic networks, but so far warm quantum memories have been susceptible to noise 
at the single photon level [3]. This problem is circumvented in cold atomic ensembles [IHH], but 
these are bulky and technically complex. Here we demonstrate controllable, broadband and efficient 
storage and retrieval of weak coherent light pulses at the single-photon level in warm atomic caesium 
vapour using the far off-resonant Raman memory scheme [TJ. The unconditional noise floor is found 
to be low enough to operate the memory in the quantum regime at room temperature. 



In an age of ever-increasing, global information trans- 
fer, there is growing demand for secure communication 
technology, such as could be provided by photonic quan- 
tum communications networks [8 . Currently, the biggest 
challenge for such networks is distance. Over short dis- 
tances, photons, interacting only weakly with their envi- 
ronment, easily and reliably carry quantum information 
without much decoherence, but intercontinental quan- 
tum communication will require quantum repeaters em- 
bedded in potentially isolated locations, because photon 
loss rises otherwise exponentially with distance [1] [2] . In 
general, these repeaters will require some sort of quan- 
tum memory, a coherent device where single photons are 
reversibly coupled into and out of an atomic system, to 
be stored, possibly processed and then redistributed. In 
order to be practically useful, this will need to have suffi- 
ciently large bandwidth, high efficiency and long storage 
time, with multimode capacity, and a low-enough noise 
level to enable operation at the quantum level. Further- 
more, to be truly practical, a critical feature of such mem- 
ories will be ease of operation near room temperature for 
use in isolated, potentially unmanned locations. 

Warm atomic vapours of alkali atoms such as rubid- 
ium and caesium are potentially excellent storage media. 
At moderate temperatures the vapour pressure is suffi- 
cient to render the doublet lines optically thick, providing 
a strong atom-photon interaction in a centimeter-scale 
vapour cell. With the addition of a buffer gas to slow 
atomic diffusion, collective atomic coherences on the or- 
der of milliseconds have been observed [9]. There is no 
need for laser trapping, high vacuum or cryogenic cool- 
ing, because atomic vapour cells can be miniaturised [TU] , 
integrated onto chips using hollow-core waveguide struc- 
tures [TT] and mass-produced. Cost-effective integration 
is likely to be crucial for quantum photonics to develop 
into a mature technology. 

A number of experiments using electromagnetically in- 
duced transparency (EIT) have demonstrated efficient 
storage and recall of optical signals in warm atomic 
vapour [12], and conditional storage and recall of her- 
alded single photons [13]. In EIT, however, the frequency 



of the signal matches the atomic resonance, and collision- 
induced fluorescence at this frequency makes the uncon- 
ditional (i.e. unheralded) noise floor of this protocol too 
high for quantum applications [3J. 

These issues are avoided in cold-atom experiments |14| - 
116) . where atomic collisions are suppressed or eliminated. 
But the technical complexity of these experiments, which 
require the trapping and cooling of atomic clouds un- 
der high vacuum, makes large-scale deployment of this 
type of memory rather challenging, especially outside 
controlled laboratory conditions. Similarly, solid-state 
memories are promising candidates for efficient and low- 
noise photon storage (T7J [18] , but currently they must be 
operated at cryogenic temperatures. 

By contrast, ensemble- type Raman memories which 
are tuned far from resonance [19] provide a possible path 
to fulfilling all of the requirements for quantum-ready 
operation at room temperature. The far off-resonant 
Raman interaction results in: (i) extremely broadband 
capability, allowing to interface the memory with con- 
ventional parametric downconversion sources |20j : (ii) 
the ability to optically switch the memory in and out 
of the quantum channel, or alternatively set the stor- 
age level to 50%, providing a straightforward nonintcr- 
fcrometric path to creating light-atom entanglement; (iii) 
very weak fluorescence noise which is predominantly non- 
synchronous with the short signal pulse; and (iv) memory 
efficiency which is insensitive to inhomogeneous broad- 
ening, allowing room-temperature operation and a path 
towards integrated implementations. Apart from the ad- 
vantage of simplicity [21] , operating at room temperature 
also makes it easier to achieve larger optical depths, and 
hence higher efficiencies []33 H2] . 

Recently, we implemented a Raman memory in a 
hot caesium-vapour atomic ensemble, demonstrating ex- 
tremely broadband, coherent operation under room- 
temperature conditions [TJ. In this letter, we address 
the remaining requirements for a memory to be able 
to function usefully in a genuine quantum application. 
Specifically, we demonstrate total memory efficiencies of 
> 30%, and with only moderate magnetic shielding, we 
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FIG. 1: (a) A-level scheme for Raman memory. A is the de- 
tuning from the atomic resonance, (b) Experimental setup. 
At ti the single photon level signal is mapped by a strong 
write pulse into a spin-wave excitation in the caesium vapour 
cell. At £2 a strong read pulse reconverts the excitation into a 
photonic mode. After polarization filtering (Pol) and spectral 
filtering with Fabry-Perot etalons, the retrieved signal is de- 
tected with single Photon counting modules (SPCM). Vertical 
polarization is indicated as (£) and horizontal polarization as 
(*+)• 



show storage times of up to 4 /is, around 2,500 times 
longer than the duration of the pulses themselves. This is 
already sufficient to improve heralded multiphoton rates 
from parametric down-conversion sources |23| . Finally, 
we make a detailed analysis of the unconditional noise 
floor of the memory, which is found to be < 0.25 pho- 
tons per pulse; that is — low enough for quantum ap- 
plications. The combination of these results shows that 
we have a quantum-ready memory, capable of handling 
quantum information in a simple room-temperature de- 
sign. 

The heart of the quantum memory is a caesium- vapour 
atomic ensemble, prepared in a vapour cell heated to 
62.5°C, that makes use of the 852-nm D2 line, with the 
6 2 <Si/2 hyperfine states serving as the ground |1) and 
storage |3) states (Fig. [l|a)). At t 1 , the signal gets 
mapped by the strong write field into a collective atomic 
spinwave excitation, and at t2, a strong read pulse re- 
converts the atomic coherence back into the photonic 
mode (Fig.jljb)). The signal is separated from the strong 
write and read fields via spectral and polarisation filter- 
ing, which is particularly important when operating at 
the single-photon level, and then detected using a silicon 
avalanche photodiode (see Methods for details). 

In the current experiment, we demonstrate total ef- 
ficiencies around 30% and a memory lifetime around 
1.5 /xs (Fig. [2]), which is double the efficiency and more 
than two orders of magnitude improvement in lifetime 
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FIG. 2: Memory efficiency. The lifetime of the memory is 
about 1.5 /is. The blue dots indicate experimental data mea- 
sured with a fast APD. Error bars represent the standard 
error in the mean. The solid line is the theoretical dephasing 
predicted for a constant magnetic field of 0.13 ± 0.05 Gauss, 
which can be attributed to the residual of the Earth's mag- 
netic field (Supplementary Information). Note that this de- 
phasing could be compensated using spin-echo techniques, or 
suppressed with improved magnetic shielding. 



over the values obtained in our previous experiment [7]- 
These results were obtained with only moderate magnetic 
shielding and residual static fields are the main dephasing 
mechanism (see Supplementary Information). In princi- 
ple, magnetic dephasing can be eliminated by improved 
shielding, enabling storage times limited by atomic dif- 
fusion to several hundred /zs [24]. However, because this 
memory has such a broad bandwidth, this measured life- 
time already corresponds to a time-bandwidth product 
of ^2500. Indeed, with memory efficiencies of 20% at 
1 /xs, this memory could already be used to improve her- 
alded photon-pair rates with typical parametric down- 
conversion sources with heralded single photon rates of 
~1 MHz. As discussed in Ref. [7], the memory efficiency 
is restricted mainly by control field power and the less ef- 
ficient, but experimentally simpler, forward-retrieval con- 
figuration [19]. 

Next, we measure the memory output using a single- 
photon counting module (SPCM) to test its performance 
at the single-photon level. Figure |3][a) shows storage at 
t = ns and retrieval at t = 780 ns later for an average 
input signal of 1.6 photons per pulse (< 1 retrieved pho- 
ton per pulse), while Figs[3jb,c) show the storage and re- 
trieval processes in detail. The strong reduction in trans- 
mitted signal (Fig. |3jb)) between the control field being 
off (green) and on (blue), and the significant amount 
of retrieved signal (Fig. |3jc), blue) demonstrates that 
the quantum memory operates well at the single-photon 
level. However, to properly characterise the noise in these 
signals, it is critical to measure the unconditional noise 
floor: the signal detected when memory retrieval is trig- 



3 



2000 - 



1 1000 ■ 





(a) 80 








60 








40 








20 













10 20 30 40 50 60 70 80 






600 



NIK) 




7X1 



(e) 











50 50 100 
Pump power (mW) 



100 



50 50 100 
Pump power (mW) 



FIG. 3: Single-photon-level data, (a) Incident (ti ~ ns) 
and retrieved (<2 ~ 780 ns) pulse. The inset is a zoom of 
the retrieved pulse showing a fluorescence noise tail. (b,c) 
Zooms around t\ and fe, showing the: incident signal (trans- 
mitted/retrieved signal with no control field) (green); trans- 
mitted/retrieved signal with the control field (blue); and noise 
(control field only) (red). Histograms are accumulated over 
360,000 runs, (d) Blue shaded area: Optical pumping on the 
3) <H> 1 2} (blue) transition. Red shaded area: Optical pump- 
ing on the |1) «-> 1 2} (red) transition. The red points are data 
and the black line is a theoretical prediction of noise due to 
spontaneous Raman scattering. The red solid line shows the 
predicted noise due only to Stokes scattering (e) Predicted 
fractional noise contributions from Stokes (black curve) and 
anti-Stokes (red curve) light. Theoretical details can be found 
in the supplementary information. 



gered without any stored input signal (Figs J3jb,c), red 
curves). In our case, this is the detected signal when the 
control field is sent in with no input signal, and is mea- 
sured to be 0.25 photons per pulse. Already substantially 
less than a single photon, this enables operation in the 
quantum regime, even with 30% memory efficiency 

We now explore the origin of the noise photons to de- 
termine if the observed noise could be further reduced. 
The inset to Fig.^a) shows that the retrieved signal peak 
is followed by an exponential tail, which is fluorescence 
noise from the excited state, |2): this process is excited 
by the control field, even far from resonance, in the pres- 
ence of atomic collisions. The red curve is a least-squares 
fit to the tail of the experimental data and yields an es- 
timated lifetime of 32±2 ns (expected value 30.5 ns [25] : 
error bars derived from multiple experiments). While 
such collision-induced fluorescence limits the usefulness 
of other more narrowband room-temperature memories 



(as shown in [3]), in our extremely broadband memory, 
the 30 ns time scale of these emissions is much longer 
than the duration of the readout event (set by the 300 ps 
pulse width and detector timing jitter). By time gating 
the detection, fluorescence therefore made a negligible 
contribution to the measured, instantaneous noise floor 
(Figs§b,c)). 

To understand where the instantaneous noise comes 
from, we then investigated its dependence on optical 
pumping. The blue shading in Fig.|3jd) indicates optical 
pumping on the |3) «-> |2) ("blue") transition, whereas the 
red shading represents optical pumping on the |1) f-> |2) 
("red") transition. Increasing the pump power on the 
blue transition partially suppresses the noise, although it 
rapidly levels off at an average of 0.25 photons per pulse, 
while the noise level rises linearly with increasing pump 
power on the red transition. 

These observations are well-described by a simple noise 
model (supplementary information) based on sponta- 
neous Raman scattering (SRS) PHUT], in which Stokes 
and anti-Stokes photons scattered during the control 
pulse are both transmitted through our etalon-based fil- 
ters (sec Methods section) and detected as noise. Since 
the model does not include collision-induced fluorescence 
[3] or leakage of the control, this suggests that these 
contributions are negligible. Furthermore, as shown in 
Fig. [3^d) and (e) , around 60% of the noise affecting the 
quantum memory is emitted at the anti-Stokes frequency 
(when operated with maximal blue pumping) and could 
be removed using further spectral filtering. This would 
already bring the true unconditional noise floor down to 
0.1 photons per pulse resulting in an unconditional signal 
to noise ratio (SNR) of 10:1 for single photon retrieval. 
The remaining signal-frequency (Stokes) noise remains 
even if the optical pumping and spectral filtering are 
perfect and has its origin in four-wave mixing seeded by 
spontaneous anti-Stokes scattering. However, even this 
noise can be eliminated if the anti-Stokes channel can 
be suppressed or rendered much weaker than the Stokes 
channel, for instance by operating the memory closer to 
resonance, so that the anti-Stokes detuning is relatively 
much larger. Finally, although our fibre-coupled detec- 
tion system is optimised for the mode of the signal, the 
noise is scattered mostly into the control mode [551 HH] j 
which suggests that the noise floor can be further re- 
duced by angle tuning the control field [30] ■ These re- 
sults show that the Raman memory, therefore, represents 
a genuine quantum-ready optical memory, functioning at 
room temperature. 

In conclusion, we have demonstrated a coherent broad- 
band single photon level optical memory at room tem- 
perature with a detailed investigation of the uncondi- 
tional single-photon noise floor. We have also extended 
the memory storage time far enough that it could be 
used to improve heralded multiphoton rates from para- 
metric down-conversion sources — an important step to- 
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wards making this type of memory useful for quantum 
repeaters. This work shows that the far off-resonant Ra- 
man memory scheme makes it possible to implement a 
quantum ready memory in warm atomic vapour, and to- 
gether with technological advances in microfabrication 
of vapour cells and hollow-core waveguide structures em- 
bedded on chips, opens a path to an integrated, scalable 
quantum information technology architecture. 



Author contribution 

KFR and PM carried out the experimental work and 
took the data. KFR, JN and NKL were responsible for 
the experimental design and theoretical analysis, and 
KCL and IAW contributed to the experiment. The 
manuscript was written by KFR with input from JN, 
NKL and IAW. 



Methods 



The atomic ensemble is initially prepared in the ground 
state |1) by optically pumping with an external cavity 
diode laser (ECDL). For the storage process, the signal 
and write pulses, both derived from a ThSapph laser os- 
cillator and an EOM [7] , are sent temporally and spatially 
overlapped into the vapour cell (Fig. Ilfb)). The pulses 
have a duration of 300 ps corresponding to a bandwidth 
of 1.5 GHz. The bandwidth of the memory is defined 
dynamically by the strong write field which is about 10 9 
times brighter than the single photon level signal field. 
The necessary filtering is provided by a Glan-laser po- 
lariser and three air-spaced Fabry-Perot etalons with a 
free spectral range of 18.4 GHz and transmission window 
of 1.5 GHz, giving a total extinction ratio of 10~ n . Note 
that the free spectral range is twice the 9.2 GHz Stokes 
shift, so both Stokes (signal) and anti-Stokes (noise) fre- 
quencies pass through the filters. As discussed in the 
main text, anti-Stokes noise could be eliminated using 
more selective filters. The detection system consists of 
single photon counting modules (SPCM) combined with 
a time-to-amplitude converter (TAC) and a multichannel 
analyser (MCA), which allows photo-detection events to 
be time-binned with a resolution of ~ 125 ps. All exper- 
imental single photon level data are histograms derived 
from 360,000 experiments. The memory operates at a 
repetition rate of 3 kHz and the storage time is set to 
800 ns. 
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SRS noise model 



The dominant noise process in our Raman memory ap- 
pears to be spontaneous Stokes and anti-Stokes scatter- 
ing during the strong control pulse. To model the noise, 
we consider the Maxwell-Bloch equations describing the 
collincar propagation of Stokes and anti-Stokes field am- 
plitudes As, Aas through a A-type ensemble with reso- 
nant optical depth d along the z-axis, in the presence of 
a control pulse, whose shape is described by the time de- 
pendent Rabi frequency f2 = f2(r), where r is the time in 
a reference frame moving with the control. In the adia- 
batic limit that the bandwidth and Rabi frequency of the 
control are much smaller than the detuning of the fields 
from resonance, the propagation equations take the form 



d, 



r s 

r 



A, = - 



AS 
1 1 



^AS _ 



1 AS 



B. 



B. 



AS 



B = -V^yO* 



Pi 

Fs 



P3_ 

1 s 

r 



(i) 



-4 a 



/it 

^AS' 



AS. 



Acknowledgements 

This work was supported by the EPSRC through the 
QIP IRC (GR/S82716/01) and project EP/C51933/01. 
KFR was supported by the Marie-Curie-Network 
EM ALL PM was supported by FASTQUAST. IAW was 
supported in part by the European Commission under 
the Integrated Project Qubit Applications (QAP) funded 
by the 1ST directorate as Contract Number 015848, and 
the Royal Society. 



where B is the amplitude of the spin wave, Fs,as = 
7 — iAs,AS is the complex detuning of the Stokes (anti- 
Stokes) field, 7 is the homogeneous linewidth of the 
|1) -H> 1 2) transition, Aas — As + S, the Stokes detuning 
Ag is equal to the detuning A of the signal field in Fig.[T] 
and 5 = 9.2 GHz is the Stokes shift. The z-coordinate 
has been normalised so that it runs from to 1. We 
have defined p\ (p^) as the fraction of atoms initially in 
state |1) (|3)). Since the noise process is weak, we assume 
that these populations, as well as the control pulse, are 
unaffected by the interaction. 
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In general, the solution is given by 



^AS,out(T) 



/oo 
^s(r,r')^s,in(r')dr' 
-OO 
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+ / L s (t, z)B in (z)dz, 
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— OO 
1 



+ / Las (t, z)Bl(z)dz, 
Jo 

where the subscripts 'in' ('out') describe the amplitudes 
at the start (end) of the interaction, and where the inte- 
gral kernels -Ks,AS) Gs,as an d is, as are Green's functions 
that propagate the input fields to compute the output 
fields. Both the Stokes and anti-Stokes frequencies are 
passed by our spectral filters, so our noise signal is calcu- 
lated by adding the average number of photons scattered 
into both Stokes and anti-Stokes modes, 

S = I (4,out(^S,out(T))dT 

(4 S ,out(^AS,out(T))dT. (2) 

oo 

It can be shown that the averaged photon numbers do 
not depend on the shape of the control pulse, but only 
on its energy, through the quantity W = |^(t)| 2 dr 
[T§] . as is usual in the transient regime of sponta- 
neous Raman scattering [35]. The initial state used 
to evaluate the expectation values is the vacuum state, 
with no Stokes/anti-Stokes photons and no spin wave 
excitations. Using the bosonic commutation relations 
[-4A,ASan(T),^A,AS,m( T ')] = S(t - t'), and noting that 
B oc |1)(3|, so that (Bl n (z)B in (z')) = p 3 5(z - z') and 
{Bi n {z)B- m {z')) — p\5{z — z'), we obtain 

/OO c oo 
/ dr'dr{|Gs(r,r')| 2 + |G AS (r,r')| 2 } (3) 
-CO J —CO 

/OO pi 
/ dzdr {p 3 \Ls (t, z)\ 2 + Pl \L AS (t, z)\ 2 }. 
-oo J a 

This expression depends only on the Green's functions, 
and can be evaluated by solving the system Eqs. ([T]) nu- 
merically, treating the amplitudes As, Aas, and B as 
classical complex-valued functions. The variation of the 
populations p±_ 3 with optical pumping power P is mod- 
elled by setting 



MP) = 
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P/Ps 



PUP, 



where P s is the saturation power and with the convention 
that negative values of P indicate pumping on the blue 
1 2) «-» 1 3) transition. The observed spontaneously scat- 
tered Raman noise is then given by S ^ S (P) = kS(P), 
which incorporates a scaling parameter into the simple 
one-dimensional theory to account for the fact that the 
Raman scattering process has a broad spatial distribu- 
tion [28] [29] , whereas the signal is emitted into a single 
narrower spatial mode. The parameter k defines the over- 
lap between the scattered mode and the detected mode, 
which is defined and filtered by a single-mode fibre opti- 
mised on the signal mode. For our experiment we have 
d = 1900, 7 = 16 MHz, A = 15 GHz and W = 30 GHz. 
To correctly describe our noise measurements, we find 
P s = 84 mW and k = 0.12 using a least-squares fit. 



Magnetic dephasing 



We model the magnetic dephasing by considering the 
evolution of the Raman coherence excited in the memory 
under the influence of a static magnetic field. Suppose 
that an atom is initially prepared in the Zeeman state 
\Fi,m{) with probability p m . within the initial hyperfine 
manifold with Fi — 4. Storing a signal pulse creates a 
spin wave, represented by applying the operator S' to 
the initial atomic state, where S — al + /3S, with / 
the identity operator, a, j3 two coefficients quantifying 
the amplitude of the spin wave, whose values are not 
important, and E the transition operator given by 



E = 



Ft 



Fi 

E ^ 

mi— — m.f — — Ff 



G(m;,TOf) |-Fi,mi)(.Ff,mf| 



(5) 



Pl (P)=l-p 3 (P), (4) 



Here C(m- U m{) is the coupling coefficient between the 
initial state and the final Zeeman state |Ff,mf), with 
Ff = 3. To model the dephasing, the absolute coupling 
strengths are not important; the relative strengths are 
computed in a straightforward manner from the Clebsch- 
Gordan coefficients describing the allowed transitions, 
once the polarizations of the control and signal fields 
are fixed. Over a time t, the atomic spins undergo Lar- 
mour precession due to the magnetic field. If the field 
has strength B in a direction parameterized by the polar 
and azimuthal angles 9, </>, the precession is described by 
the operator U = W ER, where R rotates the quantiza- 
tion axis to align it with the field, and E accounts for the 
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accumulation of phase by each Zeeman level, 



R = E ^ 

Ft 



sin <b—X\ cos 



|Fi,mi)(Fi,mi| 



+ e l{Y!Sin *- X!COS W\F t ,m t ){F h m{\\ 

mt=-F{ 
Fi 

E = J2 e im ^* Bt miXFi.mil 

mi- — Fi 
Fi 

+ J2 e imf9 ^ BS *|F f ,m f )(F f ,m f |. 
mt=-F{ 

Here X^f , Y^f are the x and y components of the spin an- 
gular momentum operators for the spins F^f, and g- lt f are 
the ^-factors determining the size of the Zeeman split- 
ting within each hyperfine manifold. In our case we have 
9i = 1/4 and = — 1/4. [ab is the Bohr magneton. 

The Raman coherence acts as a source for the re- 
trieved signal field in the presence of the retrieval control 
pulse, so the retrieval efficiency is proportional to |(S)| 2 , 
where the expectation value is evaluated just before re- 
trieval. Incoherently summing the contributions arising 
from atoms starting in each of the initial Zeeman lev- 
els, we compute the retrieval efficiency according to the 
formula 



ry(i)c< Prn^Fum^SU^US^F,, 



(6) 



--F; 



The only non- vanishing terms in the above expression are 
all proportional to |a/3*| 2 , confirming that the absolute 
values of these coefficients simply determines the scaling 
of m not its shape, so that we obtain 



Fi 

E 

mi=-Fi 



p mi \(F umi \U^U^ \F umi )\ 



The plot in Fig. [2] is produced by plotting rj(t), multiplied 
by an appropriate scaling factor, assuming a uniform ini- 
tial thermal distribution p m = l/(2Fi + 1), where the 
magnetic field B ss 0.13 Gauss and orientation (0 = 30° 
from the vertical, i.e. from the direction of the con- 
trol field polarization, and 4> — 25° from the direction of 
propagation of the optical beams) are determined from a 
least-squares fit to the experimental data. 
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